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ABSTRACT 

We present results on the spectral decomposition of 118 Spitzer Infrared Spectrograph (IRS) spectra from 
local active galactic nuclei (AGN) using a large set of SpitzerflRS spectra as templates. The templates are 
themselves IRS spectra from extreme cases where a single physical component (stellar, interstellar, or AGN) 
completely dominates the integrated mid-infrared emission. We show that a linear combination of one template 
for each physical component reproduces the observed IRS spectra of AGN hosts with unprecedented fidelity 
for a template fitting method, with no need to model extinction separately. We use full probability distribu¬ 
tion functions to estimate expectation values and uncertainties for observables, and find that the decomposition 
results are robust against degeneracies. Furthermore, we compare the AGN spectra derived from the spec¬ 
tral decomposition with sub-arcsecond resolution nuclear photometry and spectroscopy from ground-based 
observations. We find that the AGN component derived from the decomposition closely matches the nuclear 
spectrum, with a 1-cr dispersion of 0.12 dex in luminosity and typical uncertainties of ~0.19 in the spectral 
index and ~0.1 in the silicate strength. We conclude that the emission from the host galaxy can be reliably 
removed from the IRS spectra of AGN. This allows for unbiased studies of the AGN emission in intermediate 
and high redshift galaxies -currently inaccesible to ground-based observations- with archival Spitzer /IRS data 
and in the future with the Mid-InfraRed Instrument of the James Webb Space Telescope. The decomposition 
code and templates are available at http: //www. denebola. org/ahc/deblendIRS 
Subject headings: galaxies: active - infrared: galaxies - infrared:spectroscopy 


1. INTRODUCTION 

In the 5.5 years duration of the Spitzer (tWerner et al.l 
20041) cryogenic mission, the Infrared Spectrograph (IRS; 
[Houck et al.ll2004l) obtained mid-infrared (MIR) spectra with 
unprecedented sensitivity for more than one thousand AGN 
dLebouteiller et al.l 1201 it) , including quasars, Seyferts, LIN- 
ERs, radio-galaxies, blazars, etc. However, the limited spatial 
resolution of IRS observations (3.6” and 10.5” for the Short- 
Low (SL) and Long-Low (LL) modules, respectively) implies 
that the extended emission from the host galaxy usually con¬ 
tributes a significant fraction of the flux within the IRS aper¬ 
ture. Because of the difficulty in disentangling the emission 
from the AGN and its host galaxy, the shape of the AGN spec¬ 
trum at MIR wavelenghts usually has large uncertainties, ex¬ 
cept for some quasar-luminosity AGN ( Mullanev et al.ll201 It 
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Shan g et al.ll20TltLShi et al.lfloTl . 

Sub-arcsecond resolution observations with ground-based 
MIR spectrographs on 8-m and 10-m class telescopes probe 
nuclear regions in scales <100 pc for local AGN, effectively 
resolving away most of the emission from the host (see 
Figure QJ. This is evidenced by the very small equivalent 
widths (EW) of the PAH features in the nuclear spectra 
(iRoche e t al. 11200 61: iHonig et al .1120 1 (1 iGonzalez-Martfn et al.l 


20U| Escjuej ii etal. 20141 iAIonsoTferrero et al.l 120141 

Ramos Almeida et al.ll2014l) and the tight correlation between 


the MIR luminosity of the nuclear point source and the 
bolometric luminosi ty of the AGN deriv e d from hard X-ray 
observations (e.g. 

iGandhi et al.l 12009; 


Krabbe et al l 12001 iHorst et al.1 [2008; 
Asmus et al.l 1201 4R However, ground- 


based MIR observations are limited to specific windows 
of atmospheric transmission and suffer from much lower 
sensitivity. Current samples of AGN with ground-based MIR 
observations are limited to a few dozens of local sources with 
~8-l 3i im sp ectroscopy (e.g. iGonzalez-Martm et al.l 12013b 
Esquei e t al.l 12014b and a few hundred s with N-band and/or 
Q-band photometry (lAsmus et al.ll20l4h . 

An alternative approach is to separate the spectral 
components corresponding to emission from the AGN 
and the host in the integrated spectrum. Multiple 
spectral decomposi tion methods have been appl i ed to 


Lutz et al. 2004; Saiina et al. 

2007 

Hernan-Caballero et al. 

2009; Alonso-Herrero et al.l 

2012al), broadband spectral 


Mull anev et al.ll20Tl : iFeltre et al.l l2013). 


2008t iMateosetail 20151). and combi nations of both (e.g. 


Libraries of theoretical models abound in the literature both 
for the IR emission of AGN (e.g. iPier & Kro 1 i kl[T992[ H99l 
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Granato & Danese 1994c Efstathiou & Rowan-Robinson 

1995; Fritz et all :2006; 

Nenkova et al. 

2008) and galaxies 

(e.g. Silva etal. 1998 

Popescu et al 


2000; Gordon et alJ 

2001; Siebenmorgen & Krligel 2007 

Draine & Li 200'S 


suitable for fitting spectral energy distributions spanning a 
wide wavelength range, in general they cannot reproduce 
accurately, for example, the diverse shapes of the 1 Oitm 
silicate feature or the relative strength of individual PAH 
bands. More detailed models of the MIR spectrum con¬ 
taining dozens of physically motivated co mpon ents, such as 
those produc ed by CAFE ^Marshall et al.ll200~7l) or PAHFIT 
(ISmith et alJ 120071) and its variants (e.g. iGallimore et alJ 
120101) . allow to reproduce these features with high fidelity. 
However, the large number of free parameters involved (more 
than 100 in the case of PAHFIT) implies that solutions are 
often degenerate except in very high signal-to-noise ratio 
(SNR) spectra. 

Another approach that requires few free parameters yet 
provides acceptable fits for many MIR spectra is spectral 
decomposition with empirical templates. Most previous at¬ 
tempts at decomposing MIR spectra with empirical templates 
have relied on two or three sp ectral components r epresented 
by un i versal templat e s (e.g . iLaurent et alJ I2000t iTran et alJ 
200 It iN ardini et al.l l2008t IHernan-Cabal lero etal. 2009t 


Hernan-Caballero & Hatzi minaogloul 1201 It IVeilleux et al.l 
20091; iRuiz et alJ 120131) . This technique usually models 
the effects of extinction by assuming a dust screen with a 
particular extinction law that applies to at least the AGN 
component, and uses the value of Ay as an extra free pa¬ 
rameter. This simple approach has a caveat; while the MIR 
emission from the host galaxy is relatively uniform across 
galaxies, it is simply not possible to reproduce the full range 
of SEDs observed in AGN-dominated spectra by just varying 
the amount of extinction applied to a single AGN template. 
This is not criti cal when only a narrow spectral ran ge is 
considered (e.g. INardini et alJl2008t IRuiz et afl 120131) . but 
it quickly becomes an issue as the fitted spectr al range 
grows. For example, IHernan-Caballero et al.l ( 2009 ) showed 
that AGN with a steeper 5-15 pm slope compared to the 
chosen AGN template require increasingly high contributions 
from a starburst continuum (HII region) template to obtain 
acceptable fits, even if no signs of star formation (as in PAH 
emission) are observed. 

A complete description of the spectral components in a 
given spectrum requires not only finding the parameters of 
the best fitting model, but also their uncertainties. Uncer¬ 
tainties are caused by both noise in the data and degenera¬ 
cies in the template set. While most template fitting rou¬ 
tines account for the former, the latter are usually ignored. 
In the last few years, algorithms based on Bayesian inference 
have been developed to find the marginalised probability dis¬ 
tribution of physical parameters from comparison of the ob¬ 
served data with all the models in the library. Some remark¬ 
able examples are MAGPHYS ( Ida Cunha et al.l :2008l) and 
BayesCLUMPY dAsensio Ramos & Ramos Almeidal 120091) . 
However, all these methods rely on grids of theoretical mod¬ 
els with a uniform sampling of the parameter space to fa¬ 
cilitate the calculation of probability distributions. Recently, 
iNoll et aD (120091) proposed a new Bayesian-like marginalisa¬ 
tion method that is suitable for models with an irregular cov¬ 
erage of the parameter space, laying the groundwork for a 
complete description of spectral components based on purely 
empirical templates. 



Figure 1. Sketch illustrating the influence of spatial resolution in the spectra 
of local AGN. The red box represents the aperture size of a typical ground- 
based nuclear spectrum, while the green one represents the much larger aper¬ 
ture of the Spitzer/IRS Short-Low module. The spectrum observed by IRS 
is the sum of the ground-based spectrum and the spectrum in the area within 
the red and green squares, which only contains emission from the host. 

In this paper we take a new approach to spectral decompo¬ 
sition that uses libraries of observed high SNR spectra as tem¬ 
plates. We show that the method produces excellent fits that 
accurately reproduce most features in high SNR spectra, yet it 
is robust against degeneracies even in poor SNR spectra. We 
obtain reliable expectation values and uncertainties for AGN 
properties by computing full probability distribution functions 
(PDFs) using a marginalisation method that overcomes the 
inhomogeneity in the template set. By comparing the de¬ 
rived AGN properties with those obtained from ground-based 
high spatial resolution (0.4”-0.8”) MIR photometry and spec¬ 
troscopy, which exclude circumnuclear emission, we demon¬ 
strate that this kind of spectral decomposition is accurate at 
separating the emission from the AGN and its host galaxy and 
recovering the shape of the AGN spectrum. In §2 we present 
the method of spectral decomposition and our approach to the 
calculation of the expectation values of observables. §3 de¬ 
scribes the set of templates and §4 the selection of the test 
sample. In §5 we present the results of the spectral decompo¬ 
sition and perform a careful validation by comparing the AGN 
fluxes, spectral indices, and silicate strengths with those ob¬ 
tained from ground-based MIR photometry and spectroscopy. 
§6 discuses the strengths and limitations of the method in light 
of the presented results, and §7 summarizes our main conclu¬ 
sions. Throughout this paper we use a cosmology with Ho = 
70 km s -1 Mpc~', Qm = 0.3, and Qa = 0.7. 

2. THE METHOD 

Our approach to spectral decomposition relies on the as¬ 
sumption that the spectral shape of the AGN and host galaxy 
emissions that contribute to the integrated spectrum can also 
be found in other sources where either the AGN or the host 
completely dominates the emission. The emission from the 
host is itself a combination of that from the stars and the in¬ 
terstellar medium, with the former dominating in passively 
evolving galaxies (at least at A< 15 pm) and the latter in star¬ 
forming ones. Accordingly, any IRS spectrum could be re¬ 
produced by a linear combination of three IRS spectra se¬ 
lected from sets containing only extreme cases of ‘pure- 
AGN’, ‘pure-stellar’, and ‘pure-interstellar’ spectra. For clar¬ 
ity, we refer to these single-component spectra as templates. 

This procedure could raise several concerns. One is 
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that AGN feedback could have an impact on the interstel¬ 
lar medium, modifying the spectrum of the host, e.g., by 
suppressing the PAH features. However, nuclear spectra 
of nearby AGN hosts at high spatial resolutions show no 
clear sign of PAH suppression, even in the close vicinity 
of the AG N (lEsauei et~al~l[2014^ lAlonso-Herrero et al.ll2014t 

iRamos Almeida et al.l]2014fclGarcfa-Bernete et al.ll2015l) . An¬ 

other issue could be that the host galaxy modifies the intrin¬ 
sic AGN spectrum via foreground extinction from dust in the 
line of sight to the AGN. While foreground extinction can in¬ 
deed be very strong in luminous^ infrared galaxies_£LIRGs), 
mergers, and edge-on galaxies (i Goulding & Alexanderi[2009i 

lAlonso-Herrero et al.ll201 lUGonzalez-Martm et al.l2Q 13j) , we 

do not intend to recover the intrinsic spectrum of the AGN, 
but the one that would be observed at higher spatial resolu¬ 
tion. Since the extinction caused by the host is independent 
on its angular size or the luminosity of the AGN, it should be 
possible to find a ‘pure-AGN’ template with the same intrin¬ 
sic spectrum and the same level of (host induced) extinction 
but negligible host emission. 

Modelling foreground absorption in the MIR with an ex¬ 
tinction law and a few empirical or semi-empirical ‘in¬ 
trinsic’ AGN templates provides great flexibility. How¬ 
ever, this relies on the assumption of the extinction law 
and is a path already well explored (e.g. [Tran et all 1200 it 
Lutz et all l2004t IS noon et all l2004t iHernan-Caballero et alJ 
2009: iMullanev et all 1201 ilk Our aim in this work is to 
produce spectral decompositions that are completely model- 
independent. Therefore, we require the template set to be 
large and varied enough to reproduce both the diversity in the 
intrinsic spectra of AGN and hosts and the different levels of 
foreground extinction to be found. 

Our decomposition method is implemented in the De- 
blendlRS routintQ written in the IDL language. We build a 
library of templates (see © containing IRS spectra of galax¬ 
ies whose emission is strongly dominated by the stellar pop¬ 
ulation (stellar templates), the interstellar medium (PAH tem¬ 
plates), or the active galactic nucleus (AGN templates). All 
the spectra (the ones to be decomposed and the templates) are 
de-redshifted and resampled to a common wavelength grid 
A„. We choose to resample with a wavelength resolution 
A/l=0.1p,m. While the results of the decomposition do not 
depend significantly on A/l, this value is larger than the reso¬ 
lution of the original spectra and therefore increases the SNR 
per resolution element. The flux uncertainties cr(A) are calcu¬ 
lated for the resampled spectra assuming that random errors 
for adjacent pixels are not correlated. 

For every galaxy in the sample we try spectral decomposi¬ 
tions using every possible combination of one stellar template, 
one PAH template, and one AGN template. We represent each 
of these combinations with the tern of indices (i,j,k) for the i- 
th stellar template, the /'-th PAH template, and the k-th AGN 
template. The decomposition model for the spectrum F(/l) 
with the tern (i,j,k) is then: 


fjA'O - aff ar (A) + bff AH (A ) + cf AGN (A ) (1) 

where f [ s,ar (A), f^ AH (A), and f AGN (A) are the restframe flux 
densities of the templates, and the coefficients a, b, c are ob¬ 
tained by minimising the y 2 '- 

11 Available at http://www.denebola.org/ahc/deblendIRS 


2 y (F(A n ) - fij, k (A n )) 2 

x h °^ 


( 2 ) 


Since the model is linear on the coefficients, we can find 
the values that minimise x 2 algebraically. This is much faster 
than an iterative fit and allows to exhaustively explore all the 
possible template combinations in a few seconds. If any of 
the fit coefficients turns out to be negative, we fix it at zero 
and repeat the calculation for the other two until we arrive at 
a solution where all the coefficients are non-negative. 

We compare the minimum^ 2 for each combination of tem¬ 
plates to determine the one that produces the best x 2 of all 
models. We note that while the absolute minimum ofx 2 iden¬ 
tifies the best fitting model, its value is not a reliable indicator 
of the similarity between the model and the actual spectrum of 
the sour ce, since it also depends on the SNR of the spectrum 
(see e.g. [Hernan-Caballerol201 A . At high SNR, small differ¬ 
ences between the spectrum and the model can be statistically 
significant given the small photometric errors, and it is not 
unusual to obtain a x 2 per degree of freedom^ 2 = y^ 2 /(N-3) 
» 1. However, in lower SNR spectra differences between the 
spectrum and the model are blurred by the noise, and ;y 2 < 1 is 
typical. 

A complementary statistic on the differences between the 
model and the spectrum is the coefficient of variation of the 
root-mean-square error (CVrmse): 


CVrmse = 


Vstl {F(A„) - fiJ,kUn)f/N 


Z* 1 F(An)/N 


(3) 


which represents the ratio of the typical residual to the mean 
flux. This is formally equivalent to a xl that assumes flux 
uncertainties are proportional to fluxes, and therefore ignores 
the SNR of the spectra. Accordingly, a comparison of the two 
statistics provides a reliable diagnostics of the quality of the 
fit irrespective of the SNR of the spectrum. 

At constant SNR, both^ 2 and CVrmse decrease for increas¬ 
ingly accurate fits. However, unlike^- 2 , CVrmse increases if 
the SNR is degraded. This is because larger flux uncertain¬ 
ties increase the residuals (the numerator in Eq. 3), but not 
the average flux. Obtaining high values for both CVrmse and 
xl implies that no model can correctly reproduce the spectra, 
while low xl and CVrmse indicates a good fit with small rel¬ 
ative residuals in a high SNR spectrum. However, a low xl 
with high CVrmse implies the model is compatible with the 
data just because the spectrum is noisy. While xl values are 
not directly comparable among sources due to differences in 
the SNR of their spectra, the high SNR of all the templates 
ensures that for a given source we can arrange the models by 
the quality of their fit usings 2 alone. 

The empirical nature of the templates implies that we can¬ 
not expect a homogeneous and complete sampling of the pa¬ 
rameter space defined by the physical properties of interest 
(such as the spectral index of the AGN component, a, and 
the strength of its silicate feature, S s n). Instead, the most 
common values are well represented in the library by many 
templates that are similar to each other, while extreme val¬ 
ues are underrepresented. Since the probability distribution 
function (PDF) of a physical property is usually computed as 
the sum of the probabilities of all the models in given bins 
of the parameter space, the inhomogeneous sampling acts as 
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a prior that favours the most densely sampled regions of the 
parameter space. Thefore, to obtain reliable PDFs we need 
a method that is robust against variations in the density of 
templates throughout the parameter space, which would oth¬ 
erwise favour values that are overrepresented in t he library. 
One such method is the ‘max’ method described in lNoll et alJ 
d2009h . Very briefly, we take the entire range of possible val¬ 
ues for the parameter of interest, x , and split it into same-sized 
bins. For each bin, i , we find all the models that produce a 
value for the parameter within its limits, and take the lowest 
X 1 , Xr The probability of the parameter x to take a value 
inside the bin i is then: 


-xf/2 


Pi = 


£e-V /2 


(4) 


where the normalisation ensures that the sum of the probabil¬ 
ities for all the bins is 1. Then, the expectation value for x 
is: 


(x) = ^ PP 

and its standard deviation is: 




• <*» 2 


(5) 


( 6 ) 


Obtaining a realistic PDF with the ‘max’ method requires 
the bins to be large enough for each of them to contain at 
least one model. This is guaranteed for properties that depend 
on the coefficients of the three spectral components, such as 
the fractional contribution of the AGN to the restframe 12urn 
luminosity (Lp GN ), given the large number of possible tem¬ 
plate combinations. However, for properties that depend on 
the templates themselves, such as a and .S’ s ji, obtaining a cor¬ 
rect sampling depends on the library. 

We restrict the modelling of the IRS spectrum to the rest- 
frame 5-15pm spectral range. This is because the most 
prominent PAH features (which are essential to disentangle 
AGN and host emission) reside in this range, and at A> 15 pm 
discriminating between the continuum emission of the AGN 
and the host becomes increasingly difficult. In addition, the 
fraction of the total continuum emission arising from the AGN 
decreases at longer wavelengths, due to the combined effects 
of a flat AGN spectrum and a steeply rising continuum from 
dust heated by star formation (see e.g. lMullanev et alll201 11) . 
In extended sources, this effect is exacerbated by the the larger 
aperture of the Long-Low (LL; 14-38pm; 10.5”) IRS mod¬ 
ule compared to Short-Low (SL; 5.2-14pm; 3.6”) which in¬ 
creases the contamination from the host in the latter. 

Setting the long wavelength limit at 15 pm restframe also 
implies that we can use as AGN templates sources at redshifts 
up to z=1.5 while mantaining full spectral coverage. This is 
important because the diversity in AGN spectral shapes re¬ 
quires the library of ‘pure-AGN’ spectra to be as large as pos¬ 
sible for the best results (see 50. 

3. SPECTRAL TEMPLATES 

All the templates employed here are IRS low-resolution 
spectra from CASSIS of individual galaxies with full spectral 
coverage in the restframe range 5.3-15.8pm, reliable spectro¬ 
scopic redshift, and high SNR. All the templates have me¬ 
dian SNR per pixel >10 in the original spectra and >15 after 
the resampling. Typical values are much higher, with 75% of 
templates at SNR/pixel>100 after resampling. 


We select templates of three types, each of them repre¬ 
senting a distinct physical component of the integrated MIR 
emission of galaxies: the stellar population (stellar templates), 
the interstellar medium (PAH templates) and the AGN (AGN 
templates). We put great care into making sure that only one 
of these components contributes significantly to the observed 
spectrum of the galaxies selected as templates. For this, we 
rely on three criteria: the slope of the MIR continuum, the 
strength of the PAH features, and the optical spectroscopic 
classification. 

For the ‘stellar’ templates, we select 19 local elliptical and 
SO galaxies. To ensure they have negligible star formation, we 
require the PAH bands to be very weak or absent, with equiv¬ 
alent widths for the 6.2pm (EW 62 ) and 1 L3pm (EW 10 ) PAH 
bands <0.02pm (see iHernan-Caballero & Ha tziminaoglou 
mu for the details on the measurement of the EW of the 
PAH bands). We also check that the IRS spectra have a blue 
stellar-like MIR continuum and the sources are not classified 
as AGN in the NASA Extragalactic Database^ (NED). 

The 56 ‘PAH’ templates are IRS spectra of normal star- 
foming and starburst galaxies at redshifts up to z=0.14. We 
make sure these sources do not have significant stellar con¬ 
tributions to their MIR spectra by requiring both high EW of 
the PAH features (EW 62 > 1.0pm and EW 113 > 1.0pm) and a 
very weak continuum at 5 pm. We also verify that they are not 
classified as AGN in NED. 

Finally, the AGN templates are IRS spectra of sources clas¬ 
sified in the optical as quasars, Seyfert galaxies, LINERs, and 
blazars. We also include a variety of optically obscured AGN 
and radiogalaxies. The library of AGN templates includes 
sources at redshifts from z=0.002 to z=1.4 and covers sev¬ 
eral orders of magnitude in bolometric luminosity. We ensure 
that the AGN templates do not contain any significant emis¬ 
sion from the host galaxy by requiring the PAH features to 
be extremely weak or absent (EW^ < 0 . 02 pm and EW 113 < 
0 . 02 pm). 

Some of the AGN templates are from (U-)LIRGs with a 
deeply obscured nucleus, as evidenced by the strong water 
ice absorption bands and a deep silicate feature, together with 
very weak or absent PAH features. In these sources, the 
observed spectrum depends mostly on the geometry of the 
dust distribution, while the nature of the power source (ei¬ 
ther AGN, starburst, or a combination of both) plays only a 
minor role dLevenson et alJl2007t llmanishi et al.ll2007l) . As a 
consequence, we cannot evaluate the actual AGN contribution 
to the MIR emission in these sources. 

We inspect all the AGN templates with the 10pm silicate 
feature in absorption and classify them as bona fide AGN 
(AGN) or obscured nuclei (OBSCURED) depending on the 
absence/presence of H 2 O absor ption bands at ~ 6 pm and/or a 
silicate strength ,S - M | <-1.5 (see T5.2l tor the definition of S s a). 
Removing the OBSCURED templates from the set of AGN 
templates would force the decomposition routine to choose 
AGN templates that provide poor fits to obscured sources. 
Therefore, we opt to keep the OBSCURED templates, but 
consider as upper limits the AGN fractions and AGN lumi¬ 
nosities obtained in cases where the best fitting model in¬ 
volves an OBSCURED template. The final selection of AGN 
templates contains 181 spectra (159 bona fide AGN and 22 
OBSCURED sources). 

Figure [2] shows all the templates selected for use in the 
spectral decomposition, separated by type. Bona fide AGN 

12 http://ned.ipac.caltech.edu 
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Figure 2. The entire set of templates normalised at 7pm restframe and grou 
spectrum for each group. 

and OBSCURED sources are shown separately for clarity. 
Spectral shapes show very small source to source variations 
among stellar templates, as is expected for spectra domi¬ 
nated by the Rayleigh-Jeans tail of stellar emission. PAH and 
OBSCURED templates show also small variations, mostly 
caused by the strength of the PAH and silicate features, re¬ 
spectively. The largest diversity of shapes is found among 
the bona fide AGN templates. These show a wide range of 
spectral slopes, silicate feature strengths, and intensity of fine 
structure lines. 

4. SAMPLE SELECTION 

To demonstrate the power of our decomposition method 
we select a sample of local active galaxies for which there 
is available Spitzer/IRS spectroscopy as well as ground-based 
MIR high angular resolution imaging. 

Our parent sample is drawn from the subarcsecond- 
resolution MIR imaging atlas of local AGN presented in 
lAsmus et al.1 1 2014 1. The atlas includes all the publicly avail¬ 
able (by early 2014) MIR imaging with ground-based 8 -m 
class telescopes. A total of 895 independent photometric mea¬ 
surements with a variety of filters in the N and Q atmospheric 
windows are presented for 253 distinct sources, with a median 
redshift of z=0.016. 

We performed a search by coordi nates in the Cornell Atlas 
of Spitzer/IRS Source{3 QCASSIS: iLebouteiller et al.ll201 ll) 
version 6 . We used a 10” search radius and found IRS spectra 
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by categories (see text for details). The thick solid line represents the average 

for 157 of the 253 sources. We verified that coordinates of the 
nuclear point-source derived from ground-based observations 
are within the IRS Short-Low aperture in all cases. These 157 
sources comprise our main working sample. Their redshift 
distribution ranges from z=0.00024 to z=0.23, with a median 
value of ;=0.017. 

lAsmus et al.l ( 2014 ) provide optical spectroscopic classifi¬ 
cations for the sources in their sample. Our subsample con¬ 
tains 46 type 1 AGN (Syl, Syl.2, Syl.5), 8 intermediate 
type AGN (Syl. 8 , Syl.9), 64 type 2 AGN (Sy2), 13 low- 
ionization nuclear emission line regions (LINER), 24 star- 
burst/AGN composites (Cp) and 2 additional sources with un¬ 
certain classifications. 

Ground based N-band spectroscopy (-8-13iim) with spec¬ 
tral resolutions between 0.4” and 0.8” is also available for 
28 of our sources from the sample of local Seyferts com¬ 
piled by Esauei et alJ (1201 41). The o bservations where origi¬ 
nally p ublished bv lRoche et all (120061 Circinus), |Ro che et all 
<120071 NGC 5506, NGC 717 21. lYonng et al l <120071 IC 
50631. iDfaz-Santos et all <f2(TToL NGC 5135, NGC 7130), 
Honi g et all (120101 ESQ 323-G77, Mrk 509, NGC 2110, 
NGC 3227, NGC 378 3 NGC 4507, NGC 72 13, NGC 7469) 
lAlonso-Herrero et al.l d20TTl NGC 4151), Colling! (1201 ll 
NGC 2 992, NGC 438 8. NGC 53471. lAlonso-Herrero et all 
(12012bl NGC 1365), iGonzalez-Martm et al.l (2013l NGC 
1386, NGC 3081, NG C 4945, NGC 5128 , NGC 5643, NGC 
7479, NGC 7582), and lSales et all (f20T3l NGC 1808, NGC 
3281). 


http://cassis.astro.comell.edu 
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Figure 3. Comparison between the coefficient of variation of the root-mean- 
square error and the reduced x 2 in the best-fitting model for each source. Stars 
represent sources with ground-based MIR photometry and spectroscopy, 
while circles represent those with just photometry. The colour coding in¬ 
dicates the median signal to noise ratio of the resampled IRS spectra. The 
dashed line marks the CVrmse=0-1 threshold used to distinguish good and 
poor fits. 


5. RESULTS 

5.1. Spectral decomposition models 

We obtain spectral decompositions for the 118 sources 
with ground based sub-arcsec resolution MIR photometry and 
Spitzer/\RS spectra that were not selected as AGN templates. 
For the remaining 39 sources, their best fitting model is triv¬ 
ially one with themselves as the AGN template and 0% PAH 
and stellar components. We confirm that the PAH and stellar 
components are actually negligible in these sources by tem¬ 
porarily removing their spectra from the template library and 
finding their new best fitting model. The ‘max’ method allows 
to obtain expectaction values and uncertainties for physical 
parameters also for these sources. 

The parameters of the best decomposition model for all the 
sources in the sample are listed in Table |T| Columns 3 to 8 
identify the PAH, AGN, and stellar templates involved in the 
best fitting model, and list the fractional contributions of each 
of them to the integrated 5-15 pm luminosity. Columns 9 
and 10 indicate the reduced^ 2 of the model, xl, and the co¬ 
efficient of variation of the root-mean-square error, CVrmse- 
For the 118 sources that are not templates themselves, values 
of^ 2 are between 0.16 and 126 (median value: 2.17). Val¬ 
ues of CVrmse range from 0.013 to 0.51 with a median of 
0.048. The latter implies that typical fitting residuals are ~5% 
of the flux density of the spectra, equivalent to the residuals 
that would be obtained with a ‘perfect’ model in spectra with 
SNR/pixel=20 after resampling. 

Figure [3] represents CVrmse versus^ 2 for the best fitting 
decomposition model of each source. The sources with the 
highest SNR in their IRS spectra (>50 after resampling) con¬ 
centrate in a diagonal sequence, while all the sources with 
SNR<50 are found above the sequence, at increasing dis¬ 
tances with decreasing SNR. Because of the sensitivity con¬ 
straints for ground-based spectroscopy, sources in this sub¬ 
sample have on average higher fluxes and SNR in their IRS 
spectra compared to the general sample, but their CVrmse and 
xl are similar to those from other sources with the same SNR. 

Visual inspection of the fits confirms that the accuracy of 
the decomposition model decreases along the diagonal se¬ 
quence, from the bottom left to the top right corners of the 
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Figure 4. Comparison of the fractional contributions of the stellar, PAH, 
and AGN components to the restframe 5-15pm luminosity in the best fit¬ 
ting decomposition model for the IRS spectrum of each source. Big stars 
represent sources in the subsample with ground-based MIR photometry and 
spectroscopy, while circles represent those with just photometry. 
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diagram. The best models, with xl< 1 and residuals of <2%, 
are barely discernible from the spectra they replicate. On the 
other hand, the sources with the highest values of xl (> 10) 
have large residuals in their fits, indicating that their high xl 
is not due to an underestimation of the photometric errors but 
caused by the lack of suitable templates. The most extreme 
case is NGC 4945, whose peculiar spectrum is dominated by a 
deeply obscured starburst component that is outside the range 
of variation of our PAH templates. There are a few sources far 
from the sequence with^ 2 <10 and CVrmse>0.1. All these 
spectra have much lower SNR compared to those in the se¬ 
quence. Their larger flux uncertainties explain their relatively 
moderate^ 2 values in spite of the large residuals. 

The relative contributions of the stellar, PAH, and AGN 
components to the luminosity in the 5-15 urn range for the 
best fitting decomposition model (respectively, tstr, /'pah, and 
tagn) are shown in Figure [4] Since by definition these quan¬ 
tities verify tstr + /"pah + /agn = 1, only two of the three are 
independent. 

Most of the sources, and in particular all the sources in the 
subsample with ground-based spectroscopy, concentrate be¬ 
tween the lines for tstr = 0 and 0.2. This reflects the fact that 
stellar emission usually contributes a very small fraction of 
the MIR luminosity, except for a few sources with low lumi¬ 
nosity AGN and little or no star formation (most of them ellip¬ 
tical galaxies) which lie close to the horizontal axis. However, 
the blue SED of the stellar templates implies that the stellar 
components contributes a much larger fraction of the emis¬ 
sion at the shortest wavelengths (~5-6p,m). This allows us to 
quantify the contribution from stellar emission, which would 
otherwise degrade the quality of the fits and the derivation of 
the AGN properties. The highest concentration of sources in 
Figure [4] is in the region with /'agn >0-7, which is expected 
given our selection criteria. 

Figure0shows the spectral decomposition of the IRS spec¬ 
tra for the 28 sources with high-resolution ground-based MIR 
spectroscopy. In most cases the correspondence between the 
IRS spectrum and the best fitting model is notable. It is partic¬ 
ularly notable in sources where the PAH and AGN templates 
contribute comparable fractions of the MIR luminosity, such 
as NGC 1365, NGC 2992, or NGC 7582. 

The SNR of the IRS spectra is in all cases high enough for 
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Figure 5. Best-fit decomposition models for the IRS spectra of the sources with available ground-based MIR spectroscopy. The black solid line represents 
the original IRS spectrum (resampled at A/l=0.1pm resolution). The dotted, dashed, and dot-dashed lines represent, respectively, the PAH, AGN, and stellar 
components of the best-fitting model, shown in yellow. The shaded area at the bottom of each plot represents the residual (spectrum - model). 
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the contribution of noise to the residuals to be negligible. In¬ 
stead, these residuals represent real discrepancies between the 
spectrum of the source and the decomposition model. While 
the continuum is usually reproduced with high accuracy, the 
PAH bands and fine structure lines often show larger resid¬ 
uals. This is because the relative strengths of the individual 
PAH bands and emission lines vary from source to source, 
and it is unlikely for any template to match all of them. In 
some cases, the best fitting model introduces an excessive 
amount of PAH component to compensate for an AGN tem¬ 
plate with too little emission at ~8 and ~ 12pm. This is evi¬ 
denced by negative residuals for most or all of the PAH bands 
(e.g. NGC 1386, NGC 5347). The worst decomposition mod¬ 
els are found in some spectra with deep silicate features (NGC 
7172, NGC 7479, NGC 4945), which obtain their best fits 
with OBSCURED templates. This is because most of the 
IRS spectra with such deep absorptions belong to compos¬ 
ite sources which show significant PAH features overlaid on 
the continuum emission, and therefore they are not eligible as 
OBSCURED templates. 


5.2. Measurement of AGN properties 

In addition to the fractions contributed by the stellar, PAH, 
and AGN components to the MIR emission, we use the spec¬ 
tral decomposition to derive other properties of the AGN 
emission useful for diagnostics. In particular, we obtain ex¬ 
pectation values and uncertainties for the monochromatic lu¬ 
minosity of the AGN at 6pm and 12pm restframe (Lg GN , 
Lp GN ), the spectral index of the AGN emission (a), and the 
silicate strength (.S' M |) of the AGN spectrum. 

We obtain the PDF for the fractional contribution of the 
AGN to the restframe 6 pm and 12 pm flux densities with the 
‘max’ method described in >J2] Then we multiply this fraction 
by the restframe 6pm or 12pm monochromatic luminosity of 
the source in the IRS spectrum (averaged in a 0.2pm wide 
interval) to obtain Lg GN and Lp ix , respectively. 

To obtain PDFs for a and S s n, we need to measure these 
properties in all the AGN templates first. We calculate a be¬ 
tween the restframe wavelengths 8.1 and 12.5pm, assuming 
a power-law continuum of the form f v oc /. This procedure 
mini mizes the impact of the 10pm silicate feature in the re¬ 
sults dHonig et al.ll2010l :lRamos Almeida et al.ll201 1 al). 

The silicate strength dLevenson et all l2007h ISnoon et all 
120071) is defined as: 


S 


sil 


= In 


FU P ) 

FcU P ) 


(7) 


where F(A p ) and Fc(A p ) are the flux densities of the spectrum 
and the underlying continuum, respectively, at the wavelength 
of the peak of the silicate feature. We identify this peak by vi¬ 
sual inspection (see for details [Hatziminaoglou et aLll2Q15 ). 
To estimate the continuum underlying the silicate feature we 
use anchor points at both sides of the silicate feature and in¬ 
terpolate with a power-law. 

As discussed in )}2] obtaining realistic PDFs for a and .S’ s ;i 
requires that the PDFs to be constructed with a bin width that 
is large enough to find at least one template in each bin. We 
use bin widths of 0.1 for both a and S s ;i. These are compara¬ 
ble to the typical uncertainties of a and S s ii measurements on 
individual templates, and provide an adequate sampling with 
no gaps within the normal range of variation of the two quan¬ 
tities (see Figure [6]). 

Figure[T]shows examples of the PDFs of four spectral prop¬ 



Figure 6. Distribution of silicate strengths (left) and spectral indices (right) 
in the library of bona-fide AGN templates. 
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Figure 7. Examples of probability distributions obtained with the ‘max’ 
method for several spectral properties of the galaxy NGC 3227. In clock¬ 
wise order, stalling for the top right plot, the represented properties are: the 
spectral index of the AGN component, the silicate strength of the AGN com¬ 
ponent, and the fractional contributions from star formation and AGN to the 
12pm luminosity within the IRS aperture. In all cases, the histograms repre¬ 
sent probability per bin. Bin sizes are 0.01 for L 12 AGN and SB fractions, 
and 0.1 for a and S& 1 . The dashed line and shaded areas represent, respec¬ 
tively, the expectation value and l-<x confidence interval. The red solid line 
represents the value derived from the best fitting decomposition model. 

erties for one galaxy in the sample, as well as the expectation 
values and confidence intervals derived from the PDFs. 

Table [2] lists the expectation values and 1-cr confidence in¬ 
tervals (percentiles 16 to 84 of the probability distribution) of 
Uf N , Fp GN , a, and S s a for all the sources in the sample. We 
note that, by construction, the maximum of the PDF corre¬ 
sponds to the value for the best fitting model. Therefore, the 
expectation value and the value for the best fitting model are 
usually very close to each other. The only exception is PDFs 
with multiple peaks or highly asymmetric. In any case, the 
result for the best fitting model is always within the 1-cr con¬ 
fidence interval derived from the PDF. 

We compare S s q and a in Figure [8] In sources with .S’ s ji>- 1 
the best fits are obtained with the bona-fide AGN templates. 
These sources define a tight sequence where those with a<- 
1.7 show increasingly deep silicates for steeper continuum 
slopes, while for a>-1.7 the silicate appears mostly in emis¬ 
sion and is much less dependent on the spectral index. This is 
in agreement with results from lHonig et al. ( 2010 ) in a smaller 
sample. Comparison with predictions from their 3D clumpy 
torus models (their Figure 12) suggests the relation between 
a and .S’ s ii is mostly driven by variations in the profile of the 
radial dust distribution, and consistent with a small disper¬ 
sion in the average number of clouds along the equatorial line 
of sight. On the other hand, sources with deep silicate fea- 
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Figure 8. Expectation values for ,S’ M | and a of the AGN component derived 
from the full probability distribution functions. Symbols as in Figure [4] 
Error bars represent the l-cr confidence intervals. Results for sources with 
fAGN <0.25 are not shown to improve readability. 

tures (.S'sii <-1.5) obtain their best fits with some of the OB¬ 
SCURED templates. They show no obvious relation between 
.S'sii and a , as expected if foreground absorption is causing or 
at least contributing to the depth of the silicate feature (see 
e. g. iGoulding et alTl2012h . In a large sample of ULIRGs, 
lHao et al.l ( 120071) also find no correlation between a and ,S’ s ji 
when the spectral index is defined between 5 and 15p.m. 

5.3. Comparison with nuclear photometry 

In this section we validate the AGN fractions and AGN lu¬ 
minosities derived from the spectral decomposition model by 
comparing synthetic photometry on the galaxy-subtracted IRS 
spectra with the nuclear photometry from high resolution MIR 
imaging compiled bv lAsmus et al.l (120141) in their Table 2. 

The sub-arcsecond resolution MIR images were obtained 
with a variety of instruments and filters in both the N and Q 
atmospheric transmission bands. The broadband flux of the 
nuclear source (/ nuc ) is measured with procedures that assume 
either a Gaussian profile with adjustable FWHM (limited to 
<1”) or an empirical PSF obtained from the calibration star 
(see lAsmus et al.l 120141 for the details on the data reduction 
and photometric measurements.). 

We select from the entire catalogue of MIR photometric 
measurements those corresponding to sources with available 
IRS spectroscopy and obtained with a filter in the N-band. 
We also require CVrmse<0.1 in the best fitting decomposition 
model (to ignore spurious results from bad fits), and the nu¬ 
clear source to be detected at the 2<x level in either the Gaus¬ 
sian or PSF photometry. There are 456 individual measure¬ 
ments from 113 distinct sources meeting these criteria. 

To obtain the synthetic photometry we convolve the IRS 
SF spectra with the corresponding filter response curve for 
each of the ground-based measurements. We perform mea¬ 
surements in both the original IRS spectrum (/ tot ) and the 
galaxy-subtracted AGN spectrum (/agn)- The latter are ob¬ 
tained by subtracting the stellar and PAH components of the 
best fitting model from the original IRS spectrum. 

lAsmus et al.l (120141) found that PSF photometry of the nu¬ 
clear source is sometimes unreliable due to variation in atmo¬ 
spheric conditions between the observations of the source and 
the calibration star. They concluded that Gaussian photometry 
is preferable except for clearly extended nuclei. 

We can use the sources that, according to the decompo- 



Figure 9. Comparison between MIR fluxes from synthetic photometry ex¬ 
tracted from the IRS spectra (/irs) and nuclear fluxes obtained extracted with 
the Gaussian profile (/ nu c) for sources with /*agn<0.95. IRS fluxes are ex¬ 
tracted for both the integrated host+AGN spectrum (open symbols) and the 
AGN component (solid symbols). Vertical dotted lines connect measure¬ 
ments corresponding to the same source and filter. The diagonal line marks 
the 1:1 relation. Big open squares and circles mark measurements for NGC 
1097 and NGC 7130, respectively (see text). The color coding indicates the 
fractional contribution from the AGN to the 5-15 pm luminosity derived from 
the spectral decomposition. Horizontal bars indicate the photometric errors 
of nuclear measurements. The inset plot shows the distribution of the ratios 
between AGN (filled histogram) or total IRS (open histogram) and Gaussian 
nuclear fluxes. 

sition, have a very strong AGN component (tagn^O- 95) to 
test the flux cross-calibration of the Spitzer/IRS and ground- 
based observations. Since in these sources the host galaxy 
contributes <5% of the flux in the IRS aperture, we expect 
/nuc ~ /agn ~ /tot- We find that Gaussian fluxes are tightly 
correlated with /agn- Their ratios have a symmetric distribu¬ 
tion with median 0.998 and l-cr dispersion 0.08 dex or 20%. 
This is roughly consistent with estimates of the uncertainty in 
the absolute flux calibration of 5-10% for the IRS spectra and 
10-15% for the ground-based photometry. Ratios calculated 
from PSF fluxes, on the other hand, show higher dispersion 
(0.17 dex) and a bias towards / nU c</AGN (median ratio: 1.22). 
We note that the distribution of ratios for PSF fluxes is also 
symmetric, suggesting that extended nuclear emission is not 
the main culprit of their higher dispersion. 

Figure [9] compares / tot and /agn from the IRS spectrum 
with Gaussian nuclear fluxes from the ground-based photom¬ 
etry for the remainder of the sample, that is, the sources with 
r agn <0-95. Dotted lines connect values of /, ot (open symbols) 
and /agn (solid symbols) corresponding to the same nuclear 
photometric measurement, and represent the decrease in the 
flux density caused by the subtraction of the stellar and PAH 
components from the IRS spectrum. We have color-coded the 
symbols according to the tagn of the sources. 

The correlation between total IRS fluxes and nuclear fluxes 
is not as tight as in the /'agn>0.95 subsample, and there are 
important deviations (always in the sense / to t>/nuc) in many 
of the sources with tagn<0.5. The distribution of / 0 t//nuc 
(empty histogram in Figure|9]l peaks at 1:1 but is highly asym¬ 
metric, with a long tail for ratios >1. This is the expected re¬ 
sult, since host emission arises mostly from outside the region 
sampled by nuclear photometry. However, when we consider 
the host-subtracted AGN spectra, the tight correlation with 
the nuclear fluxes is restored, and it holds remarkably well 
even for sources with AGN fractions as low as rAGN~0.3 (cyan 
symbols). 

The distribution of /agn/Ziuc in the i-agn<0.95 subsample 
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Figure 10. Comparison of AGN components derived from the best fitting de¬ 
composition model with nuclear photometry and spectroscopy from ground- 
based observations for the outliers in the relation shown in Figure [9] NGC 
1097 (left) and NGC 7130 (right). The solid area represents the observed IRS 
spectrum and the solid line its best fitting decomposition model. The dotted 
line represents the AGN component. Diamonds with error bars represent the 
nuclear photometry. The magenta line represents the nuclear spectrum of 
NGC 7130. 


(filled histogram in the inset plot in Figure 0 is not com¬ 
pletely symmetric, and it is slightly biased towards /agn < 
/nuc (median value: 0.895). This could indicate a system¬ 
atic underestimation of ~10% in the AGN fluxes from the 
spectral decomposition. However, it could also indicate that 
some of the nuclear emission (10% on average) arises in the 
host galaxy, not the AGN. The detection of PAH features in 
the sub-arcsecond resolution nu clear spectra of some of these 
sources favours the latter (see $5. 4I >. The boldest such case, 
marked with big open circles in Figure [9] is NGC 7130. This 
galaxy hosts a Compton thick AGN and a nuclear s tarburst 
ilLevenson et all 12005b iGonzalez-Delgado et al.l[l998h that is 
known to contribute a significant fraction of the MIR emission 
even on sub-arcsecond scales (lAlonso-Herrero et all 120061 
iDfaz-S antos et aH 120101) . Interestingly, the spectral decom¬ 
position finds that the AGN contributes ~60% of the nuclear 
luminosity in the 8-12|rm range (Figure[l0] right panel). 

Nearly all the highest values of /AGN//nuc correspond to 
measurements from a single galaxy, NGC 1097 (marked with 
big open squares in the figure). This galaxy has a low lu¬ 
minosity LINER or Seyfert 1 nucleus as well as a y oung, 
compact (<9pc) nuclear star cluster dMason et alj2007l) . The 
IRS spectrum from CASSIS also includes extended emission 
from a kpc-scale circumnuclear star formation ring. It is strik¬ 
ing that the best decomposition model fits remarkably well 
the IRS spectrum (Figure [TO] left panel), yet it predicts AGN 
fluxes 2-5 times higher than observed for the nuclear source, 
which could itself be dominated by emission from dust heated 
by the star cluster instead of the AGN dMason et all 120071) . 
This source has one of the lowest AGN fractions in the sam¬ 
ple (Vagn=0.20). At such low fractions, the relative errors in 
the predicted AGN fluxes are large. While the 1-cr confidence 
interval for the 12 pm flux density of the AGN is 60-130 mJy 
(compared to ~25 mJy for the nuclear source), the PDF is 
wide and it estimates at 5% the probability of /agn<25 mJy. 

Excluding these two outliers, the 1-cr dispersion in 
/AGN//nuc for the galaxies with rAGN<0.95 is 0.12 dex or 32%. 
Since the dispersion in the /'agn >0.95 subsample is 0.08 dex 
(20%) due to photometric uncertainties alone, we estimate a 
25% uncertainty in the fractional contribution of the AGN to 
the observed IRS fluxes. This is in good agreement with a 
mean value of 23% for the relative error in the 12pm AGN lu¬ 
minosities derived from the decomposition, and confirms that 
PDFs derived with the ‘max’ method provide reliable uncer¬ 
tainties. 


5.4. Comparison with ground-based sub-arcsecond MIR 
spectra 

There are 28 sources where we can directly compare the 
galaxy-subtracted AGN spectrum from the IRS spectroscopy 
with a sub-arcsecond resolution ground-based spectrum. 

The absolute flux calibration of the ground-based nuclear 
spectra has larger uncertainties compared to both the photom¬ 
etry from ground-based imaging and IRS spectra mostly due 
to slit losses. To overcome this, we have re-scaled the nuclear 
spectra to match the available nuclear photometry. Scaling 
factors are calculated as the ratio between fluxes from the nu¬ 
clear (Gaussian) photometry and synthetic fluxes measured 
on the nuclear spectrum with the same filters. Since most 
sources have several photometric measurements, we obtain a 
weighted average of all the scaling factors for a given source, 
with weights: w,- = l/oj, where <x ; is the uncertainty in the 
scaling factor for the i’ h measurement. 

Figure [IT] compares the AGN spectra derived from the de¬ 
composition with the nuclear spectroscopy and photometry 
from the ground. With few exceptions, there is a remarkable 
correspondence between the shape of the observed nuclear 
spectra and those predicted by the decomposition model. This 
correspondence is particularly interesting in sources with a 
strong PAH component in the IRS spectra, such as NGC 1365, 
NGC2992, NGC 3227, or NGC 7469, where it would be very 
difficult to constrain the AGN luminosity (let alone its SED) 
by other means. There are, however, significant discrepancies 
in the normalization of many sources, with nuclear fluxes up 
to ~30% lower or higher compared to IRS AGN fluxes. While 
uncertainty in the absolute calibration of both the IRS spectra 
and the nuclear photometry contribute to these discrepancies, 
it is insufficient to explain the extreme cases. In the following 
we discuss other interpretations for the origin of these dis¬ 
crepancies that apply to subgroups of sources. 

Some sources with ground-based nuclear fluxes signifi¬ 
cantly (>20%) higher than IRS AGN fluxes have clear PAH 
emission features in their nuclear spectra that evidence nu- 
clear starbursts. We l l known cases are NG C 1808 (e.g. 
iTacconi-Garman et al.1 [l996t ISales et all l20l3h . NGC 3227 
([Davies et al.ll20061l2007HHonig et alJI20m and NGC 7130 
( Dfaz-Santos et al.1 20101) In these sources, the nuclear star- 
burst contributes a fraction of the unresolved nuclear MIR 
emission, which is therefore higher than the emission from 
the AGN alone. Accordingly, AGN fluxes derived from the 
decomposition might be better estimates of their actual AGN 
output. On the other hand, if no high angular-resolution spec¬ 
trum is available, the spectral decomposition does not provide 
any information as to whether there is nuclear star formation 
or not. 

Other sources show a slight nuclear excess compared to the 
IRS AGN flux in spite of very weak or absent PAH features in 
the nuclear spectrum (NGC 1386, NGC 4507). An inspection 
of the spectral components and residuals in Figure [5] reveals 
that the AGN component derived from the decomposition is 
probably underestimated in these two sources. Their residu¬ 
als resemble an inverted PAH spectrum, indicating the PAH 
component has been boosted at the expense of the AGN to 
improve the fit of a pronounced 8 pm bump. 

In NGC 3783 the nuclear spectrum is slightly above the 
original IRS spectrum. This is unphysical, since the IRS aper¬ 
ture is much wider than the ground-based one. Inspection of 
the pointing information in the IRS spectrum reveals that the 
nuclear source was slightly off-center in the IRS slit. 
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Figure 11. Comparison of IRS decomposition results with the ground-based photometry and spectroscopy. The solid black and blue dashed lines represent the 
restframe 7.5-13.5pm section of th e IRS spectrum a nd the AGN component, respectively. The green diamonds represent the Gaussian fluxes for the nuclear 
point source in the catalogue from lAsmus et alj <20141) . with error bars representing 1-cr uncertainties. The magenta solid line represents the full-resolution 
ground-based spectrum scaled to fit the Gaussian photometry. The light magenta shaded area represents the 1-cr uncertainty in the scaling factor for the nuclear 
spectrum. 
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On the other hand, the cases where the AGN component 
in the IRS spectrum is above the nuclear spectrum could be 
explained by dust heated by the AGN on arcsecond scales 
(i.e. dusty ionisation cones). The Gaussian photometry from 
lAsmus et all d2014l) does not recover all the nuclear flux if 
the source has a FWHM>1” or the profile is less concentrated 
than a Gaussian. This seems to be the case at least for NGC 
3081, NGC 4151, and NGC 5128, where Gaussian fluxes are 
higher than PSF fluxes, yet lower than IRS AGN fluxes. 

While uncertainties in the absolute flux calibrations and the 
physical extension of the AGN-heated MIR emitting region 
introduce significant dispersion between ground-based and 
IRS fluxes, the shape of the AGN SED should be relatively 
unaffected and therefore provides a more demanding test of 
the decomposition model. 

To quantify the differences in the SED of the nuclear point 
source (as represented by the ground-based spectrum) and the 
AGN spectrum derived from the decomposition model, we 
use the spectral index of the continuum (a) and the strength 
of the 10pm silicate feature (5 s n) as indicators. 

We calculate a for the nuclear spectra usi ng the same 
method applied to the AGN templates (see ^5. 2b . with anchor 
points at 8.1 and 12.5urn. For the two sources whose nuclear 
spectra are truncated shortwards of 12.5pm, we use 8.1 and 
12.0pm instead. We estimate the uncertainty in a via error 
propagation assuming a 5% uncertainty in the fluxes. 

To estimate S s n in the nuclear spectra we place the first an¬ 
chor point for the underlying continuum at 8.1 pm and the sec¬ 
ond somewhere between 12.7 and 13pm, putting great care 
in avoiding noisy regions and the 12.81pm [Ne ii] line. As 
we did for the AGN templates, we identify the peak wave¬ 
length of the silicate feature by visual inspection. Error bars 
for S s ii are estimated by performing Monte Carlo simulations 
in which flux densities are varied within their uncertainties. 
We note that these errors account for the uncertainty in .S' s ;i 
due to noise in the spectra, but do not include the system¬ 
atic error introduced by the selection of the continuum anchor 
points and the uncertainty in the wavelength of the peak of the 
silicate feature. 

Figure [TJ] shows the spectral index of the ground-based 
nuclear spectrum (a nuc ) versus the expectation value for the 
AGN component in the IRS spectrum (o'agn). There is a 
good correlation between the two, albeit with a significant 
dispersion. The largest discrepancies correspond to sources 
with very weak AGN components (NGC 1808, NGC 7582) 
or poorly fit decomposition models (NGC 4945, NGC 7172). 
For the remaining sources a nuc and ctagn values are compati¬ 
ble within their uncertainties. The median value of the differ¬ 
ence in spectral indices, Ictagn - awl. is 0.19. 

Figure[13]compares the S s n from the nuclear spectra and the 
AGN component in the IRS spectra. There is a good agree¬ 
ment between the two, even in the sources with the deepest 
silicates (NGC 4945, NGC 7172) which obtain poor fits with 
large residuals. The median absolute difference in .S' s ;i mea¬ 
surements is 0.10, with a small bias towards smaller nuclear 
,S' s ji values compared to those from the AGN component on 
the IRS spectra. 

Obtaining such tight correlations for both a and S s u con¬ 
firms that the spectral decomposition successfully separates 
the AGN and host galaxy components of the emission, and 
that the AGN properties derived are reliable (within their un¬ 
certainties). We emphasize that even at the sub-arcsecond res¬ 
olution of ground-based spectroscopy, nuclear spectra do not 
always represent accurately the AGN emission, which can be 



Figure 12. Comparison between spectral indices measured directly on the 
nuclear spectra and those for the AGN component derived from the spectral 
decomposition of the IRS spectrum. The dashed line marks the 1:1 relation. 



AGN S sII 


Figure 13. Comparison between silicate strengths (S s n) measured directly 
on the nuclear spectra and those for the AGN template in the best fitting 
decomposition model for the IRS spectrum. 

extended or overlap with a nuclear starburst. This is, arguably, 
one of the factors contributing to the scatter in Figures [12] and 
[T3l which would be even smaller if a better proxy for the ac¬ 
tual AGN spectrum could be found. 

6. DISCUSSION 

We have found that the main limitation of our decomposi¬ 
tion method is the existence of categories of sources whose 
spectra are not well reproduced, at least with the current li¬ 
brary of templates. We have identified galaxies with deep sili¬ 
cate absorptions as one important such category. We consider 
that the reason for the poor fits in galaxies with deep absorp¬ 
tions is the small number of suitable templates in the library, 
which is caused by the difficulty in finding pure-AGN and 
pure-PAH spectra with deep absorption features within the 
limited (albeit large) set of available IRS observations, since 
many of the spectra with such deep features belong to com¬ 
posite sources. The relatively small size and restricted selec¬ 
tion criteria of our test sample implies that we cannot rule out 
other categories of sources which could also prove difficult to 
model with the current library. Even for source types that typ¬ 
ically obtain good fits, there probably exist extreme cases of 
exotic spectra whose properties are beyond the normal range 
of variation covered by the templates. However, we expect 
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these to be few in numbers. 

Degeneracy was also a legitimate concern, since the large 
number of available templates suggested there could be many 
combinations of them capable of producing good fits with 
comparable^- 2 values for any particular source. However, we 
found that while this is typically the case, the fractional con¬ 
tribution to the MIR emission from each of the three spectral 
components is usually well constrained. The reason for this 
remarkable lack of degeneracy is probably the singular char¬ 
acteristics of the 5-15 pm spectral range, where the spectral 
shapes of the stellar, interstellar (PAH), and AGN emissions 
are clearly different from each other. 

Because the spectra for the stellar and PAH components 
show very limited variation, we can recover not just the 
strength of the AGN emission but also its spectral shape. 
Decomposition models that obtain good fits for a particular 
source usually involve AGN templates very similar to each 
other, indicating that it is not easy for the stellar and PAH tem¬ 
plates to compensate for an AGN template with, for example, 
the wrong spectral index or silicate strength. 

The presence in the library of many similar (but not identi¬ 
cal) templates is not unnecessary redundancy: they help im¬ 
prove the fits by increasing the chances of reproducing fine 
details in the spectrum, like emission lines. This is important 
because in high SNR spectra these fine details dominate over 
photometric errors in the residuals and cause large x 1 values 
if the template set is not varied enough. Having multiple sim¬ 
ilar templates does not complicate the calculation of PDFs for 
AGN properties thanks to the mechanics of the ‘max’ method. 
In fact, the resolution of the PDFs is limited by the size of 
the template set, since at least one template in each bin is re¬ 
quired. More detailed PDFs allow for better estimates of the 
uncertainties in AGN properties. 

We have compared our results from the spectral decompo¬ 
sition with ground-based high spatial resolution imaging and 
spectroscopy as the ultimate test for our method. There is a 
general agreement between galaxy-subtracted AGN and nu¬ 
clear fluxes within their uncertainties, and a good correspon¬ 
dence in a and .S' M |. This has positive implications for the 
analysis with this method of the IRS spect ra of AGN at higher 
redshifts (see iHatziminaoglou et al.|[2015i) . or those to be ob¬ 
tained in the near future with the Mid-InfraRed Instrument 
(MIRI) onboard the James Webb Space Telescope (JWST). 
However, there are also discrepancies, which in some cases 
suggest a poor decomposition and in others seem caused 
by a significant contribution from the host even at the sub- 
arcsecond scales of the nuclear photometry and spectroscopy. 

The sensitivity of the fit to the shape of the AGN template 
depends on the relative contribution of the AGN component 
to the MIR emission within the IRS aperture. When it is 
low, the uncertainty in the spectral shape of the AGN is large. 
This means that AGN properties other than its luminosity are 
poorly constrained if the host galaxy strongly dominates the 
emission within the extraction aperture. This fundamental 
limitation of the decomposition method implies that high spa¬ 
tial resolution observations are necessary if the AGN does not 
dominate the integrated emission of the galaxy. Furthermore, 
even for AGN-dominated sources, ground-based high angular 
resolution observations are still required to infer the spatial 
distribution of the host emission (for example, to distinguish 
between nuclear, circumnuclear, or galaxy-wide star forma¬ 
tion) and to detect weak fine structure lines. 

7. SUMMARY AND CONCLUSIONS 


In this paper we have demonstrated a new data-driven 
method for spectral decomposition of mid-infrared spectra 
that uses as templates sets of real spectra from individual 
sources whose emission is strongly dominated by a single 
physical component. We have shown that this decomposi¬ 
tion method produces fits with an excellent accuracy for the 
Spitzer/IRS spectra of composite sources, with residuals a 
factor ~2 lower than those of a competing algorithm with a 
higher number of free parameters. While some peculiar spec¬ 
tra obtain poor fits due to their physical properties falling out¬ 
side the range of variation available from our template set, the 
combined^- 2 and CVrmse criteria allow us to identify these 
cases and set them apart for a customized analysis. 

We have applied the ‘max’ method to obtain PDFs for 
model parameters and physical properties of the AGN. Tak¬ 
ing into account the full PDF instead of just the best fitting 
model enables an adequate treatment of degeneracy in the pa¬ 
rameters derived from the decomposition models, which al¬ 
lows for realistic uncertainties. We have obtained restframe 
6 pm and 12pm luminosities, spectral indices (a), and silicate 
strengths (5 S ;i) for the AGN spectra after removal of the emis¬ 
sion from the host. We find no clear trend of either a or S s n 
with the AGN luminosity. There is a tight correlation between 
a and 5 s n that applies to all AGN except those with deep sili¬ 
cate absorption, where strong reprocessing by dust erases the 
signature of the dominant power source. 

To validate the AGN luminosities derived from the spec¬ 
tral decomposition we have compared sub-arcsecond resolu¬ 
tion N-band photometry from ground-based observations with 
synthetic photometry derived from the IRS spectra before and 
after subtraction of the host emission. We find a good agre- 
ment between nuclear fluxes and those from the AGN compo¬ 
nent in the IRS spectrum, with a typical dispersion of 0.12 dex 
in the ratio between the two. We conclude that this dispersion 
is dominated by uncertainty in the relative flux calibration be¬ 
tween the ground-based photometry and the IRS spectra. 

We have further validated the decompositions using a sub¬ 
sample of 28 sources with ground-based high angular res¬ 
olution 8-13pm spectra. We find a remarkable correspon¬ 
dence between the shape (as described by the spectral in¬ 
dex a and the silicate strength .S/i) of the observed nuclear 
spectra and the AGN spectrum derived from the decompo¬ 
sition model, even in sources whose IRS spectrum is dom¬ 
inated by emission from the host galaxy. We interpret dis¬ 
crepant cases as caused by either extended AGN emission or 
nuclear starbursts. The decomposition of IRS spectra pre¬ 
dicts a and S s a of the nuclear spectra with a typical accu¬ 
racy of 0.19 and 0.10, respectively. These small errors con¬ 
firm that our spectral decomposition is successful at separat¬ 
ing the AGN and host component of the emission and pro¬ 
vides a reliable estimate of the shape of the AGN spectrum. 
This allows for unbiased studies of the AGN emission in in¬ 
termediate and high redshift galaxies, which are inaccesible to 
current ground-based instruments, using archival Spitzer/IRS 
observations and in the near future with JWST/MIRI. 
The decomposition code and templates are available at 
http://www.denebola.org/ahc/deblendIRS 
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APPENDIX 

COMPARISON WITH DECOMPIR 

To asses the accuracy of the fits obtained with our decomposition method, we compared the residuals of the best fitting model 
produced by our method with those from another decomposition routine, DecompIR (iMullanev et all2()TTh . We chose DecompIR 
because it produces very good fits for a wide variety of Spit zer /IRS spectra and is frequently used in the recent literature (e.g. 
iGoulding et al.ll20i~2t iMagdis et al .11201 3t iDrouart et al.ll2014l) . 

For the comparison we used the Spitzer/IRS spectra of the 118 sources in our samples that were not selected as AGN templates. 
We pre-processed the spectra using the exact same procedures for de-redshifting, rebinning, and truncation of the /i> 15.8um 
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Figure 14. Comparison of the coefficient of variation of the root-mean-square error (CVrmse) for the best fitting model obtained with DecompIR and our 
decomposition routine. Symbols as in Figure[4] The dotted line represents the 1:1 relation. The inset plot represents the distribution of the ratio between CVrmse 
values obtained with the two routines (solid histogram) and its median value (dashed line). 
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Figure 15. Comparison of decomposition models for two sources with deep silicate features. The shaded area represents the rebinned IRS spectrum. The best 
fitting model for the combined host+AGN emission is shown in solid line (for DeblendIRS) and dot-dashed line (for DecompIR). The AGN spectra fitted by 
DeblendIRS and DecompIR are represented by the dotted and dashed lines, respectively. 

portion of the spectra. Because the starburst templates from DecompIR are not defined for 4 <6 pm, we also truncated this section 
of the spectra. We used the latest public version of Decomp I llQ with all nine parameters defining the model set free. We ran 
DecompIR five times for every source, each time with one of the five different starburst templates, and took for each source the 
solution with the lowest^- 2 . 

For each source, we computed the coefficient of variation of the root-mean-square error (CV RM se; see Eq. 3) for the best 
fitting decomposition model obtained with our method and DecompIR. Figure [14] compares the two values for the 118 sources. 
It shows that residuals from DecompIR are usually larger, with 90% of sources lying above the line that marks the 1:1 relation. 
The median ratio of the CVrmse is 1.95, which indicates that the typical residuals with our routine are a factor ~2 smaller than 
those obtained with DecompIR. In some extreme cases, residuals can be up to a factor ~10 smaller (sources in the top left corner 
of Figure IT4l>. While DecompIR produces smaller residuals in 10% of the sources, nearly all of these obtain very good fits with 
both methods (CVrmse<0-04). 

Interestingly, all the sources with deep silicate features that obtained poor fits (large residuals) with our routine also obtain 
equally large or even larger residuals with DecompIR. This is despite the fact that DecompIR uses two free parameters to fit the 
extinction in the AGN and the host separately. By default, DecompIR limits the range of available optical depths at 9.7pm, T9.7, 
to the interval [0,5] for both the AGN and the host. In many of the sources with deep silicates the best fitting model uses 79 7=5 
for either the AGN or the host (or both), suggesting a wider range is needed. After increasing the maximum T9.7 to 15, some of 
the sources obtain marginally lower x 1 with higher extinction values. However, the impact on the residuals is negligible. 

In Figure[l5]we show as examples the best fitting models from DeblendIRS and DecompIR for the two sources in the subsample 
with ground-based spectroscopy that obtained the worst residuals with DeblendIRS: NGC 4945 and NGC 7172. Residuals with 
DecompIR are even larger, indicating that this kind of sources are difficult to fit also for methods that use a freely varying 
extinction coefficient. 
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Table 1 

Parameters for best-fitting decomposition model 


Source name 

z 

PAH template 

*PAH 

AGN template 

TAGN 

stellar template 

r STR 

xl 

CV rmse 

ARK 120 

0.0327 

MCG -02-33-098 

0.036 

[HB89] 1435-067 

0.875 

NGC 5831 

0.089 

0.55 

0.016 

CGCG 041-020 

0.0360 

2MASX J14311978+3534184 

0.342 

ESO 103- G 035 

0.590 

NGC 4377 

0.068 

1.41 

0.034 

Circinus 

0.0014 

NGC3187 

0.046 

IRAS 12071-0444 

0.871 

IC 2006 

0.083 

13.45 

0.070 

ESO 005- G 004 

0.0062 

ESO 353- G 020 

0.268 

NGC 6300 

0.714 

NGC 5831 

0.018 

6.80 

0.060 

ESO 121-IG 028 

0.0405 

NGC3187 

0.044 

ESO 103- G 035 

0.773 

MESSIER 085 

0.183 

0.63 

0.055 

ESO 138- G 001 

0.0091 

ARP 256 NED01 

0.055 

VII Zw 244 

0.936 

NGC 5831 

0.009 

2.14 

0.034 

ESO 140- G 043 

0.0142 

NGC3310 

0.101 

PG 1302-102 

0.863 

NGC 5812 

0.037 

0.49 

0.016 

ESO 141- G 055 

0.0371 

NGC 3801 

0.023 

3C445 

0.896 

NGC 1700 

0.080 

0.23 

0.013 

ESO 198- G 024 

0.0455 

NGC3310 

0.035 

PG 1302-102 

0.929 

NGC 4660 

0.036 

0.49 

0.027 

ESO 286-IG 019 

0.0430 

MCG +08-11-002 

0.281 

NGC 4418 

0.685 

NGC 5831 

0.034 

7.50 

0.073 

ESO 297- G 018 

0.0252 

NGC 0877 

0.253 

IC 5063 

0.643 

NGC 1700 

0.104 

1.67 

0.037 


Note. — Table 1 is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 


Table 2 

AGN properties derived from decomposition 


Source name 

z 

log LfN/erg s’ 1 

1-cr interval 

log Lf 2 GN /erg s’ 1 

1 -<x interval 

a AGN 

1-cr interval 

.S'si, AGN 

1-cr interval 

ARK 120 

0.0327 

43.98 

43.80-44.09 

44.03 

43.96 - 44.07 

-1.03 

-1.45- 

-0.70 

0.31 

0.09 - 0.44 

CGCG 041-020 

0.0360 

42.75 

42.53-42.89 

43.16 

43.06-43.23 

-2.01 

-2.59- 

-1.56 

-0.89 

-1.33--0.56 

Circinus 

0.0014 

42.52 

42.44 - 42.66 

43.11 

43.10-43.13 

-2.25 

-2.42 - 

-2.14 

-1.47 

-1.65--1.40 

ESO 005- G 004 

0.0062 

41.80 

41.72-41.85 

42.06 

42.01-42.10 

-1.71 

-2.00 - 

-1.35 

-1.53 

-1.85--1.28 

ESO 121-IG 028 

0.0405 

42.46 

42.15-42.66 

42.93 

42.85-42.99 

-2.25 

-2.87- 

-1.69 

-0.67 

-1.15--0.26 

ESO 138-G 001 

0.0091 

43.22 

43.11 -43.31 

43.47 

43.45-43.49 

-1.74 

-2.00 - 

-1.52 

0.12 

-0.09 - 0.23 

ESO 140-G 043 

0.0142 

43.32 

43.24-43.40 

43.52 

43.48-43.55 

-1.70 

-2.00 - 

-1.42 

0.21 

0.01-0.30 

ESO 141-G 055 

0.0371 

43.91 

43.75-44.00 

44.03 

43.98 - 44.06 

-1.33 

-1.77- 

-1.03 

0.27 

0.05 - 0.38 

ESO 198-G 024 

0.0455 

43.36 

43.21-43.45 

43.58 

43.53-43.61 

-1.73 

-2.12- 

-1.38 

0.23 

0.01-0.36 

ESO 286-IG 019 

0.0430 

43.38 

< 43.40 

43.83 

< 43.84 

-1.80 

-1.97- 

-1.73 

-3.99 

-4.17--3.93 

ESO 297- G 018 

0.0252 

42.51 

42.37 - 42.62 

43.01 

42.95-43.06 

-2.20 

-2.61 - 

-1.91 

-0.17 

-0.41 --0.03 


Note. — Table 2 is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 







